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ABSTRACT: 

We consider the chiral expansion for the reaction vrA^ ttttN in heavy baryon chiral 
perturbation theory. To order we derive novel low-energy theorems that compare 
favorably with recent determinations of the total cross sections for vr+p vr+vr+n and 
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1. Over the last few years, new data for the reaction vrA^ ttttN in the threshold 
region and above have become available, see e.g. [Q [|] and the compilation in 
0. The interest in this reaction stems mostly from the fact that it apparently offers a 
possibility of determining the low-energy vrvr elastic scattering amplitude whose precise 
knowledge allows to test our understanding of the chiral symmetry breaking of QCD. 
However, at present no calculation based on chiral perturbation theory is available which 
links the pion production data to the tttt vrvr amplitude in a model-independent fashion. 
Consequently, all presently available determinations of the S-wave vrvr scattering lengths 
from the abovementioned data should be taken cum grano salts. In the framework of 
relativistic baryon chiral perturbation theory, Beringer |^ has considered tree diagrams. 
In that approach, however, there is no strict one-to-one correspondence between the 
loop and the small momentum expansion due to the nonvanishing nucleon mass in the 
chiral limit [0. This problem can be circumvented if one considers the nucleons in the 
non-relativistic limit. This was first used by Gasser and Leutwyler [§] (and others) and 
later formulated in terms of heavy quark effective field theory methods in ref. 0. We will 
make use here of the two-flavor formulation detailed in ref . |jlO[ . The aim of this letter is 



to show that the first two terms in the chiral expansion of the threshold amplitudes for 
ttN —y TTTcN lead to a set of low-energy theorems which indeed can be tested against the 
available data for tt^p —>■ Tr+vr+n and ir^p —>■ n^TT^n close to threshold.^ Naturally, at the 
next stage one has to consider the following terms in the chiral expansion to make contact 
with the vrvr interaction. 

2. To be specific, consider the process vr^A^ tx^tx^N , with N denoting the nucleon 
(proton or neutron) and 'a, 6, c' are isospin indices. At threshold, the transition matrix- 
element in the vr^A^ centre-of-mass frame takes the form 



ia ■ k 



Di{t^5'"' + T^S""^) + D2 T^S''^ (1) 



where k denotes the three-momentum of the incoming pion and the amplitudes Di and 
D2 will be subject to the chiral expansion as discussed below. They are related to the 
more commonly used amplitudes A2ij^^, with / the total isospin of the initial vrA^ system 
and /tttt the isospin of the two-pion system in the final state, via 

A32 = Aio = -2D, - 3D2 (2) 

which have recently been determined p|. In what follows, we will also consider the total 
cross section for the reactions n^p n^Ti^n and n^p TC^TT^n. At present, only in these 
two channels there exist accurate data in the 20... 30 MeV region above threshold. The 
data of ref. for vr+p — ir'^n^p are still too sparse and inaccurate in the threshold region 



^This result was indirectly contained in ref. but not made explicit and is much more transparent in 
the formulation used here. 
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which we are investigating. Assuming that the amphtude in the threshold region can be 
approximated by the exact threshold amplitude, the total cross section can be written in 
a compact form, 



m 



2 



^tot(s) = '-i-^X{s,m^,M^)Ts{s)\r],D^ + V2D2\' S (3) 

with m the nucleon and the pion mass, respectively and s the total centre-of-mass 
energy squared. r3(s) denotes the conventional three-body phase space and A(x, y, z) the 
Kallen-function. The r^i 2 are channel- dependent isospin factors and S" is a Bose symmetry 
factor. For vr+p Tr+vr+n and n^p n'^ii'^n we have rji = 2\/2, r]2 = 0, S = 1/2 and 
T]i = 0, i]2 = -\/2, S = 1/2, in order. In the threshold region, one can approximate to 



a high degree of accuracy the three-body phase space and flux factor |Tl| by analytic 
expressions so that 



M2J3(2 + /i)(2 + 3//) , , 

"-^^'^^ = 128vr^(l + 2,.)n/ '^^^^ + V^D^l' S {T^ - T^^ (4) 

Here, T.^ is the pion kinetic energy in the laboratory frame, T^r = (s — ml — M^^) / {2mi) — 
M^i where the subscript '1' denotes the particles in the initial state. We furthermore 
have introduced the small parameter n = M^/m ~ 1/7. This completes the necessary 
formalism. 

3. In QCD, the chiral expansion of the amplitude functions Di and D2 takes the formQ 

D = fo + fifi+ + ... (5) 

modulo logarithms. We are interested here in the first two coefficients of this expansion. 
To calculate them, we make use of heavy baryon chiral perturbation theory as detailed in 
ref.[0. The pertinent effective Lagrangian has the form 

-^^eff = ^ItN + ^ilf + (6) 

with C^^^^ given in ref. and the standard meson Lagrangian e.g. in ref.|l^. The 



(1 2) 

diagrams with insertions from jC^^ which are non- vanishing at threshold are shown in 

fig.l. Notice that the much debated next-to-leading order vrvr interaction does not appear 

(2) 

at this order in the chiral expansion. It is important to note that from only terms 
which are kinematical 1/m corrections contribute. None of the low-energy constants ci_2,3 
related to elastic vrA^ scattering (in particular to the isospin-even S-wave scattering length 
a+) ||13[ appear at order (here, q denotes a small momentum or a meson mass). One 



^Here, D stands as a generic symbol for Di^2- 
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can therefore write down low-energy theorems for £'1^2 which only involve well-known 
physical (lowest order) parameters, 



Do 



9a 
9a 



7M. 



3 + 



2m 
17M 



f ) + OiM'J 
+ 0{Ml) 



(7) 



8F3 V" ■ 2m 

with qa the axial-vector coupling constant. In what follows, we will always use the 
Goldberger-Treiman relation = gnNF^/m to calculate the numerical values of D12 
(with g.„N the strong pion-nucleon couphng constant). 
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Fig.l: Diagrams which give the contributions to D\^2 up-to-and-including order 0{Mt^\ 



The circle-cross denotes an insertion from L 



(2) 



There are potentially large contributions from diagrams with intermediate A(1232) states 
of the type Ml / ( m-A — m — 2Mt^ ) , which numerically would be of the order 10 ■ M^.^ We 
have checked that no such terms appear from diagrams involving one or two intermediate 
A resonances. Consequently, the chiral expansion is well behaved but not too rapidly 
converging. The order Mj^ corrections give approximatively 50% of the leading term. 
As we will discuss below, the calculations of Beringer in relativistic baryon chiral 
perturbation theory indicate that further 1/m suppressed kinematical corrections are 
small. To get an idea about the corrections to eqs.(0,H) we have also calculated the 
imaginary parts of the threshold amplitudes from the one-loop diagrams shown in fig. 2. 
These start to contribute at order M3 with the result 



ImDi 



Im D2 



_V3glMl 
128 7rF5 

5VSg\Ml 
64 vr Fi 



0{Ml] 



(9) 



(10) 



•^Possible large A-contributions starting at order have yet to be investigated in a systematic 
fashion together with loop effects and alike. 
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Of course, at this order there are other contributions to the real parts of -Di,2, so one 
should consider the corrections given in eqs.(|9|,p!0D as indicative. 
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Fig.2: One-loop diagrams which give a nonzero Im D\^2 at order 0{M^). 



4. Let us now turn to the numerical results. We use Ft, = 93 MeV, gj^N = 13.4, 
m = 938.27 MeV and = 139.57 MeV. This amounts to Di = 2.4 fm^ and D2 = -6.8 
fm^ or using eq.(^ 

I32 = 2AM;', Aio 



As2 = 2AM~\ Ain = 5.5 (111 



which compare favourably with the recent determinations of ref. , A32 = 2.07±0.10 
and Aio = 6.55 ± 0.16M~^. If one assumes that the imaginary parts eqs.(p|,p!OD set the 
magnitude for the order corrections of ReDi^2; then one expects Di to change very 
little and D2 by approximatively 30%. We have also calculated the cross sections for 
TT^p —>■ TT+TT+n and 7r~p — > n^n^n using eq.(0). These are shown in fig. 3 in comparison 
to the existing data. Notice that we have calculated the matrix-elements in the isospin 
limit, in fig. 3 we have shifted the resulting cross sections to account for the corresponding 
thresholds as proposed by Beringer[^]. To a high degree of accuracy, one can parametrize 
the cross sections calculated from the first two terms of the chiral expansion of D12 by 
the simple forms using eq.(^) 

= 0.225 /ib ( ^^^) = 172.4 MeV) (12) 



<o?^^"^°"(T.) = 0.442 /xb (^ YoMeV J ^^'^'^ = ^^^'^ ^^^^ ^^^^ 

as shown by the dashed lines in fig. 3. The solid lines differ very little from the ones in 
ref. indicating that higher order 1/m corrections (which are summed up in the relativis- 
tic approach) are fairly small. This also means that the approximation of using the exact 
threshold amplitude in the threshold region is a very good one for the first 30 MeV. The 
advantage of the heavy mass approach used here is the strict one-to-one correspondence 



thr" 
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between the loop and small momentum expansion. The abovementioned expectations of 
higher order corrections from ImDi 2 are indeed such that they can improve the descrip- 
tion of the data since the first/second channel allows to test D1/D2, respectively. 



TT* p -> TT* TT* n 




[GeV] 

TT" P -> TT° TT° n 




[GeV] 

Fig. 3: Total cross sections for vr+p n'^n'^n and 7T~p n^n^n in comparison to the 
data. Squares: ref.0, diamonds: ref.|^ and octagon: ref. 0. 

5. We have considered the first two coefficients of the chiral expansion for the threshold 
ttN ttttN amplitudes and derived a set of low-energy theorems, eqs.(|^J^,p|,[lO|), which 
only involve well-known physical parameters. We have also shown that the corresponding 
cross sections agree with the empirical ones close to threshold. Of course, to go further, 
one has to consider loop diagrams as well as contributions from resonance exchange (like 
e.g. the A(1232)). Ultimately, this will tell how accurately one can in fact get to the 
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elastic tttt amplitude from data on single pion production. Work along these lines is in 
progress. 
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